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We propose a new theoretical approach to the excitation spectrum of superfluid "^He. It is based 
on the assumption that, in addition to the usual Feynman density fluctuations, there exist localized 
modes which describe the short range behaviour in the liquid associated with microscopic cores 
of quantized vortices. We describe in a phenomenological way the hybridization of those two 
kinds of excitations and we compare the resulting energy spectrum with experimental data, e.g. 
the structure factor and the cross section for single quasi-particle excitations. We also predict the 
existence of another type of excitation interpreted as a vortex loop. The energy of this mode agrees 
both with critical velocity experiments and high energy neutron scattering. In addition we derive 
a relation between the condensate fraction and the roton energy and we calculate the reduction of 
the ground-state energy due to the superfluid order. 



PACS: 67.40-w,67.40.Vs,63.20.Ls,62.60+v 
I. INTRODUCTION 

In this paper we present an alternative description of 
the energy spectrum of superfluid *He. Among the stan- 
dard approaches to this problem, we mainly find either 
the effective Hamiltonian of Bogoliuhovtl, or the varia- 
tional approach proposed by Feynmana. 

The Bogoliubov description relies on the assumption 
of a weakly interacting Bose gas. Its main success was to 
obtain a linear dispersion for the long wavelength excita- 
tions starting from the quadratic spectrum of free bosons. 
This linear phonon branch cannot be obtained from sim- 
ple perturbation theory but corresponds to a RPA-like 
description of this problem. This linear behavior is in- 
deed one of the main features of the experimental spec- 
trum (fig.l). Nonetheless this success, neither the initial 
slope of the dispersion £(fc), nor the rest of the spectrum 
beyond the linear part, is quantitatively obtained from 
this approach. _ 

A alternative description was proposed by FeynmanEI, 
based on a variational form of the wavefunction for the 
low-lying excited states which leads to an excitation spec- 
trum of the form 
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where m is the mass of the atoms and S{k) is the static 
structure factor of the liquid. |-This expression, using the 
experimental structure factom, is shown in fig.l. It re- 
produces quantitatively the low momentum linear behav- 
ior, since in this limit the f-sum rule is exhausted by 
these density fluctuationsQ. Unlike the weakly interact- 
ing Bose gas limit, it also reproduces qualitatively the 
roton minimum. Despite these successes, this variational 
wavefunction, although it includes Bose symmetry, does 
not contain special features characterizing the superfluid 
phase. For instance, since it contains only information 
about the static structure factor, it gives essentially the 



same energy spectrum (|l|) both below T\ in the superfluid 
phase and in the normal phase (T > T\) where exper- 
imental results exhibit very different behaviors for the 
spectrum of these two phases. Unlike S{k) which almost 
does not change, the dynamical structure factor S{k,uj) 
is very sensitive to the superfluid transition. This shows 
upj-in the vanishing of the roton and maxon peaks above 

In order to get a better quantitative, agreement in the 
roton region, subsequent refinementsQ'D of the original 
variational wavefunction were proposed, which include 
more localized structure on top of the delocalized density 
fluctuations (phonons). This was achieved by introduc- 
ing more variational parameters, at the expenses of the 
initial interpretation of the excitations as density fluc- 
tuations. Despite sizeable improvements for the energy 
spectrum and other measurable quantities including the 
scattering intensity, this approach does not fully answer 
the question of the nature of the superfluid behavior. 

An alternative approach to the superfluididty of ''He 
using path integral Monte Carlo calculationscl was pro- 
posed. Although it does not provide explicit results for 
the spectrum, it showed unambiguosly how superfluidity 
is related to the coherent exchange of atoms along rings 
of all sizes. 

Our motivation in order to build a phenomenological 
model for the superfluid is based on the features dis- 
cussed previously. We assume that the main features 
of superfluidity show up in the spectrum in the large 
momentum range above the linear phonon part. This 
includes the roton minimum (at energy Aland the satu- 
ration for very large fc's at the energy 2Ad. We propose 
to view superfluid "^He as the condensation in one state 
of a set of localized states with long-range interactions. 
This is in contrast to the model of hard-core bosons on a 
lattice without interactions. We shall describe these lo- 
cal modes as two- level systems (TLS) of bare energy E^). 
This approach is analogous to the description of dipolar 



glasseglil where such a picture, based on two-level tun- 
neling states, has been successfully used to describe the 
low temperature excitationsapeasured by specific heat or 
heat transport experimentsE£l. 

The TLS we consider are localized on atomic scales 
and may correspond to quantum coherent rings of a few 
atoms only, which are supposed to be responsible for the 
superfluid behaviour of the liquid. These short rings cor- 
respond to vortex cores of sizes R ~ lA, and may be 
considered as local dipoles. Large ring structures arise 
from the long-range interaction between the TLS. 

The effective Hamiltonian Hioc we propose for these 
interacting local modes can be diagonalized by a Bogoli- 
ubov transformation. The resulting spectrum bears sim- 
ilarity with the weakly interacting Bose gas expression, 
although it does not assume a weak interaction. Indeed, 
in order to recover a gapless phonon branch at low mo- 
mentum (k ~ 0), we need to consider strong interac- 
tions, namely interactions of the order of the bare energy 
£^0- Unlike ordered structures of dipoles on a lattice for 
which it is possible to calculate the interaction term, for 
liquid ^He (and for glasses), the interaction term is ob- 
tained from the experimental energy spectrum. Despite 
this shortcoming, this model gives closed analytical ex- 
pressions for quantities such as the scattering intensity of 
neutrons by single quasiparticle excitations, the relative 
variation of the ground-state energy and the condensate 
fraction in terms of a single quantity, namely the strength 
of the interaction, so that it becomes possible to relate 
them. 

In order to compare this approach to the Feynman vari- 
ational scheme, we use an equivalent formulation of our 
problem based on a description of the interaction between 
the localized modes as due to a virtual exchange of Feytu 
man phonons,— This is very much in the way HopfieldEHI 
and Andersonll3 considered the exciton problem in dielec- 
tric media. This description is equivalent to the previous 
one based on the effective Hamiltonian Hioc provided we 
use the dipolar approximation which corresponds to large 
momentum. For low momentum, both the density fluctu- 
ations and the large scale exchange cycles correspond to 
the same degrees of freedom and cannot be disentangled 
anymore. In the large momentum limit the energy spec- 
trum results from the hybridization of these two kinds 
of excitations, which represent independent degrees of 
freedom of the superfluid. This way we gain a better 
understanding of the roton part of the spectrum. In ad- 
dition we find a new localized excitation branch at energy 
E2 — SE'oj which we interpret as microscopic vortex-loop 
excitations of the superfluid. These type of excitations 
are usually thought as being completely independent of 
the phonon-roton excitations, while here we propose a 
model that describes them both together. 

This paper is organized as follows. We begin by intro- 
ducing the effective Hamiltonian of the strong coupling 
description in section II. In section HI, wc obtain expres- 
sions for the scattering intensity of the single-excitation 
branch and compare it to the experimental data. In sec- 



tion IV, we evaluate the reduction in the ground-state 
energy of the superfluid and the condensate fraction. In 
section V, we consider the dipolar approximation and 
the resulting hybridization scheme. We compare the re- 
sulting spectrum and the static structure factor to the 
experimental data. Finally, in section VI, we discuss the 
vortex-loop branch and we conclude in section VII. 



II. THE EFFECTIVE HAMILTONIAN 
DESCRIPTION 

We consider superfluid ''He as a set of local states de- 
scribed as two- level systems of bare energy Eq. Their 
Hamiltonian is then 



Kc = Y.^obk^'"^ 



(2) 



where translational invariance allows us to write opera- 
tors in /c-space. The operators bk obey bosonic commuta- 
tion relations. Thiis holds in the limit of a low density of 
localized modes lij, i.e. for the case of no multiple occupa- 
tion of a site. This condition corresponds to vortex-cores 
or small exchange rings which cannot be multiply excited 
on the same site. The value of Eq will be determined later 
on. 

The interaction between these localized-modes can 
generally be written ast3 
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The diagonal part in this Hamiltonian corresponds to 
the hopping of a local mode between sites, and the off- 
diagonal part describes the creation (or the anihilation) 
of pairs of local modes on distinct sites. The matrix el- 
ement X{k) depends on the microscopic details of the 
interaction. Fo|^-a lattice structure like excitons in a 
crystal or solidEJ bcc ^He, it is possible to express X{k) 
as a summation of dipolar terms. This is not possible in 
the superfluid phase, for which we have no lattice struc- 
ture. 

The Hamiltonian Hioc is diagonalized by the Bogoli- 
ubov transformation j3k — u{k)bk + v{k)h^ ^k- The re- 
sulting spectrum is 



E{k) = ^EoiEo + 2Xik)) 
and the two functions u{k) and v{k) are given by 
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The corresppading ground-state wavefunction for the lo- 
cal modes isEJ 
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where \vac) is the state with no local modes, i.e. anihi- 
lated by bk, namely bk \vac) — 0. The ground-state is a 
coherent state of pairs of localized-modes. This coher- 
ent state defines a global phase which, once fixed, breaks 
the gauge-invariance. The excitation given by P^k |^o) 
corresponds to the breaking of a pair of momentum fc, 
leaving an unpaired local-mode. 

The expressions of the bare energy Eq and X{k) stem 
from the following remarks. The equality E{k) = Eq is 
achieved in the zero coupling limit {X{k) — 0) and it cor- 
responds to a normal phase for the local modes charac- 
terized hy Vk — in the relations (pi) and (p). This corr&j 
sponds to the termination point of the energy spectrumEl 
for which, at large momentum fc, the energy goes to 2A 
where A is the roton energy (fig.l). Beyond the termina- 
tion point, there is no quasi-particle excitation and only a 
structureless free-recoil scattering characteristic of a nor- 
mal liquid. We shaU therefore take for the bare energy 
the value Eq = 2^. 

In the opposite limit, i.e. for low momentum fc — > 0, 
we demand that the spectrum (^) will be gapless, namely 
E{0) = 0. This requires 



X(0) = f = -A 



(7) 



Notice that for dipoles on a latticdiJ, this condition 
is a self-consistent definition of Eq since the energy re- 
quired to excite locally a dipole from the ground-state 
corresponds to the one needed to flip it with respect to 
the coherent background. This excitation changes the 
sign of the interaction energy and costs 2 | X{k — 0) \ = 
2 A = Eq, which corresponds to the previous definition 
of the energy of the local mode. 

For arbitrary k, X{k) is obtained by comparing Eq.(|j) 
with the measured energy spectrum (fig. 2). This shape 
oi X{k) is similar to the expression calculated for dipples 
on a simple cubic lattice along various directionsEJo. 
This could be helpful towards the development of a mi- 
croscopic model for the superfiuid phase. 

It is useful at this stage to compare our approach with 
the weakly interacting Bose gas limit. There, the Hamil- 
tonian is given within the Bogoliubov approximation by 

HwiBG =^^{^k'^ NaVk) a\ak 

k 

+ Y^ NaVk (4a-fc + a^a-fe) (8) 



where N^ = |(ao)| , £*: = ^ A;^/2m, and Vk is the effec- 
tive potential between the bosons at wave vector k. The 
operators a\ and ak correspond to the creation and ani- 
hilation of an atom and not of a localized- mode as in (O). 



The Bogoliubov excitation spectrum which corresponds 
to d) is 



v/£fc [ek + 2iVo^fe) 



(9) 



This spectrum is linear in the fc — > limit { E c^ hkc), 
with a sound velocity defined by c = y^NoVo/m. 

The physics described by the two Hamiltonians Hi^c 
and HwiBG is very different. The weakly interacting 
Bose gas has a non-interacting limit given by the ideal 
Bose gas which, at T=0, is known to be fully condensed 
in the state k = and which has already a broken gauge 
symmetry. The interaction depletes the condensate and 
changes the spectrum to linear at low momentum. On the 
other hand, the Hamiltonian Hioc corresponds to a set of 
independent localized modes with no broken symmetry 
in the non-interacting limit {X{k) = 0). To restore in 
the linear momentum regime both the condensation and 
broken gauge symmetry, we need a non-zero interaction 
X{k). Since this interaction is large (^(0) = —Eo/2), it 
is not a small perturbation to the non-interacting case. 
In that sense, it is a strong-coupling description of the 
superfiuid. 

We would like to point out that the entire excitation 
spectrum given by (H) is unique to the superfiuid phase, 
as only there we can define the function X{k) and local 
mode energy Eq. In particular, the phonon mode is dif- 
ferent frota the zero-sound mode that appears in the nor- 
mal fluidcj. We shall now use this approach to obtain an- 
alytical results that we shall compare subsequently with 
experimental measurements. 



III. SCATTERING INTENSITY 

The neutron scattering intensity is a direct probe of the 
density fluctuations in the liquid and may be described 
using the dynamic structure factor S{k,(jj). It is usuafe 
accepted since the work of Miller, Pines and Nozieres Ell 
that we can split the total contribution to S{k,Lu) into 
two parts. 



S{k,u;) = NZ{k)6{huj - e{k)) + S'^^\k,Lu) 



(10) 



where the first term accounts for single quasi-particle 
excitations of weight Z{k), while the second describes 
multiparticle excitations. This separation is justified 
at low temperature and low momentum, typically for 

k < O.bA . In this regime, integrating (lin) over the 
energy and noticing that S'^^\k,uj) vanishes in the low 
momentum limit we obtain Z{k) ~ S{k), which results 
also from the Feynmaii. theory. The comparison with 
the experimental data U shows that it works indeed in 
the low momentum regime mentionned above, but fails 
to describe the non linear part of the spectrum, except 
perhaps for the position of the maximum. 

Since the excitations of the local modes from the 
ground state (n) involve breaking pairs of coherent local 



modes, the probability of a neutron to scatter inelasti- 
cally on such an excitation is proportional to the occupa- 
tion density of these pairs at each momentum given by 

(^fe^^Lfe) = UkVk, so that 



Z{k) =^ Airk'^ loUkVk 



fll) 



where /q is an arbitrary normalization constant and 
where the factor Ank"^ comes from the three dimensional 
phase space. Using the expression (|^) for Uk and Vk, we 
obtain UkVk — \ eIJ) '^^li^h together with (|j) gives 



Z{k) = 7rfc2/o 



E{k) 



m-^ 



(12) 



Using the experimental E{k) in expression (|12|) we ob- 
tain a differential cross-jSection which agrees well with 
the experimental resultsu) obtained at saturated vapor 
pressure (S.V.P.) as shown in fig. 3. In the low momen- 
tum limit, we recover the proportionality between S{k) 
and Z{k) and therefore Z(k) is linear with k. We empha- 
size that as a result of both ( [l2|) and the saturation of 
E{k) to 2A at large k, the scattering intensity Z{k) van- 
ishes identically at this point. This indeed corresponds 
to the experimental measurement of Z{k), and could not 
be obtained neither from weakly interacting Bose gas de- 
scription nor by variational approaches. 



IV. DESCRIPTION OF THE CONDENSATE 

A. Condensate fraction 

In the microscopic description of a Bose liquid, the 
condensate is usually characterized by the condensate 
fraction, n^/n, which measures the relative occupation 
by the atoms of the lowest energy state. At zero tem- 
perature it is equal to one for an ideal Bose gas, while 
a finite interaction depletes the lowest energy state with 
a corresponding decrease of no/n. It is possible to cal- 
culate this fraction in terms of the interaction potential 
only within the Bogoliubov description of a weakly inter- 
acting Bose gas. More generally, it is knownE3 that this 
ratio rio/n is related to the zero momentum limit of the 
non-condensate particle distribution in the ground state, 
Uk, namely 
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where c is the sound velocity and m is the mass of a ^He 
atom. For the real liquid, we do not known how to ob- 
tain no/n in terms of other measutable quantities. This 
ratio has been obtained indirectlyo from a measurement 
of the distribution of the non-condensed-atoms. Numer- 
ically, using path integral Monte-CarlcQ similar values 
have been obtained. An analytical expression has been 



proposed which is based on the assumption that the de- 
pletion of the, condensate is related to the thermally ex- 
cited rotonscil. 

In the present model, it is also possible to obtain the 
divergent part of the ground state occupation number 
Hk — v{ky of the local modes. Using (^, it is given by 



"■fe^o = 



A 
2?ifcc 



(14) 



Since both ([T^ ) and ( [Til) describe the same behaviour, 
namely the condensation of the ^He atoms, then, by 
equating the two corresponding residues, we obtain the 
relation 



no 
n 



(15) 



This relation predicts a direct and unexpected link be- 
tween the condensate fraction no/n and the characteris- 
tics of the energy spectrum at T = 0, namely the sound 
velocity c and the roton energy A. The behavior of both 
no/n and A/{mc^) measured as a function of the pres- 
sure are in agreement with the relation (jl^) up to a factor 
of nearly two which we shall discuss later on. 

Another relevant feature of (|l^) is that it allows to re- 
late the quantity no/n which is critical at the superfluid 
transition and therefore is proportional to the amplitude 
of the complex order parameter, to the roton energy A 
which, in our model, is the modulus of the interaction 
matrix element at k — namely |X(A; = 0)| = Z^ A 
critical behavior of A was indeed recently observedEa in 
high resolution neutron scattering experiments, and sup- 
ported by Raman scattering experiment^. This seems 
to indicate that the roton energy A (and therefore the 
interaction matrix element X(k)) is related to the super- 
fluid order as it appears in (|l5|). 



B. Condensation energy 

Another quantity which characterizes the condensate 
is the condensation energy AEq- It is obtained by rewrit- 
ing Hioc in (H) in terms of the operators /3fc and Pi, 
namely 

Hioc - E Eik)plPk +J2l iE{k) ~ [Eo + X{k))) 

k k 

(16) 

The second term in ( [iq ) is the shift in the ground state 
energy AEq- Using Eq.(Q) we obtain 



^^G = -^E(^w-2Af 
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k^ {E{k) - 2Ay dk (18) 



where kmax is the largest value of the momentum and 
corresponds to the termination point of the quasi-particle 
spectrum. Vsp ^fif4^ specific volume per atom. Using 
the experimcntaBuH (fig.l) energy spectrum for E{k), 
we obtain 



AEg^-5.0±1.0K{SVP) 
AEg ci -3.0 ± 1.0K{P ^ 24:atm) 



(19) 



for saturated vapor pressure and P=24atm respectively. 
In fig.4 we plot the integrand P {E{k) - 2A) in ^. 
It shows that the main contribution to the reduction in 
the ground state energy comes from the roton part of the 
spectrum due to the larger phase space volume, compared 
with the phonon region (fc — > 0). |— ■ 

Both from experimental datau3 and path integral 
Monte-Carlo calculationsB, it is possible to deduce the 
values for the change in the kinetic zero point energy 
of the atoms between the normal phase at Tx and the 
superfluid phase in the limit T — j 



AEun^2.b±i.{)K{SVP) 
AEmu ^ 1.3 ± 3.0K{P = 24atm) 



(20) 



for saturated vapor pressure and P=24atm respectively. 
First, we notice that the large experimental uncertainty 
makes any quantitative comparison difficult. Moreover, 
the two quantities AEq and AEkin do not measure 
exactly the same property. While AEq measures the 
change in the ground state energy of the local modes 
between the non-interacting and interacting limits, both 
at T=0, AEkin measures the change in kinetic energy 
per atom between two different temperatures. Since the 
variation of tha-uolume of the superfluid between T\ and 
T = is smallcZI, it is usually accepted to atribute the 
change in the kinetic energy mainly to the condensation 
phenomenon. In our model this condensation energy is 
calculated per local mode, nevertheless we find that both 
quantities decrease for an increasing pressure. 



C. Effective mass 



Although the relations (Iq) and (18) for the conden- 
sate fraction and the condensation energy give both the 
right order of magnitude and the expected bahaviour as 
a function of pressure, they are consistently larger than 
the experimental values by a factor of about 2. This fac- 
tor may stem from the fact that our expressions are given 
as a function of the local modes while the experimental 
data is obtained per '^He atom, and there is no direct 
one-to-one correspondence between these two. There- 
fore, the mass term which appears in (^ instead of be- 
ing the bare mass m of a ^He atom is the effective mass 
rUeff of the bare local mode. From the above compari- 
son between no/n and AEq with experiments we obtain 
that rUeff ~ 2m. Although such a relation can be ob- 
tained only from a microscopic calculation beyond our 



phenomenological model, it is interesting to notice that 
a similar result was obtained numerically using path inte- 
gral Monte-Carlo calculationsQ. There, the effective mass 
m* of a tagged *He atom which does not participate in 
the Bose permutations was found to be m* ~ 2m. Since 
such a tagged atom is a local node or defect in the su- 
perfluid, it corresponds to the bare local mode we con- 
sider in the Hamiltonian Hioc- As we saw in (|5|) and (|g), 
the bare local mode of energy Eq appears for the non- 
interacting case {X{k) = and v{k) = 0). Then, it does 
not contribute to the coherent (superfluid) ground state, 
i.e. behaves as a local node of the order parameter. This 
correspondence may justify taking m,eff — 2m, bringing 
both ( p^ ) and (|lj) into agreement with the experimental 
data represented in fig. 5, where the experimental valuesE^ 
for no/n, A and c are obtained independentlynZI. 



V. HYBRIDIZATION OF LOCALIZED AND 
DENSITY MODES 

In the previous sections we described the superfluid as 
resulting from the condensation of a set of localized and 
interacting modes in the strong coupling limit. This de- 
scription, as we emphasized, is very different from the 
Bogoliubov weak coupling limit. Since it is aimed to de- 
scribe the full energy spectrum, we would like to relate 
it to the Feynman variational scheme. In this section we 
shall show that in addition to the phonon-roton branch of 
excitations around the roton momentum, there is another 
branch of excitations. We shall discuss the properties of 
this excitation in the-ilight of experiments on thermal nu- 
cleation o£-vorticefca and on large momentum neutron 
scatteringLJ. To that purpose, we propose in this sec- 
tion a description equivalent to the previous interacting 
local mode problem, where the interaction between these 
modes results from the exchange of virtual phpaons. This 
is very much in the same spirit of HopfieldEj and An- 
derson considering the exciton problem in dielectric me- 
dia. There localized excitons, taken as dipoles, interact 
through the exchange of photons. For the superfluid, the 
local modes play the role of the excitons, and interact 
through the exchange of Feynman phonons. This picture 
is only valid in the range of large enough momentum be- 
yond the linear part around /c — > 0, which corresponds 
to density excitations (i.e. Feynman phonons). It is only 
in this regime that the two types of excitations repre- 
sent independent degrees of freedom of the system. The 
Hamiltonian describing Feynman density excitations is 



iJo = ^ e{k)a\ak 



(21) 



where afe|-|are Bose operators and e{k) is the Feynman 
spectrumB given by the relation (nh, expressed in terms 
of the static structure factor S{k). The Hamiltonian de- 
scribing the coupling between these phonons and the lo- 
calized modes is 123 



Hr 



k 

h.c. 



(A(fc, Eo)bk + /i(fc, EQ)ak) (flfc^ + a^k) (22) 



and the total Hamiltonian is 



H = Ho 



Hio, 



Hr 



(23) 



Since H and Hioc provide two equivalent descriptions of 
the same problem, the functions A and /i are related 
to X{k). If we now assume that the localized modes 
are locaL jdjpolar excitations, then there is a specific 

tioiEJ-H 



relati 



between these functions, given by A(fc) 



^^0 [-^) ' and ^l{k) = -Eo^. As we shall dis- 
cuss later on, this dipolar approximation is valid only 
for large momentum, i.e. in the roton region. The two 
modes Uk and bk are independent degrees of freedom, so 
that their commutator is zero. The total Hamiltonian M, 
can be diagonalized using the canonical transformationO 

afc = A{k)ak + B{k)bk + C(fc)at_fe + D{k)b^ ^k (24) 
5fe = B{k)ak + A{k)bk + D{k)a)^k + C(fc)fe^-fc 

which describe respectively the lower and upper 
branch of the energy spectrum. The functitms 

A{k), B{k), C{k), D{k) can be written down explicitl}it3. 
We point out that the roles of the phonons and localized- 
modes are interchanged between the two branches, as can 
be seen by writing the ground-states of both branches: 



ioi)=n 
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The corresponding dispersion relation is 

[E^{k) - e^{k)] [E^{k) - E^] = -6X{k)EoE^{k) (26) 

We first notice that for a vanishing coupling X{k) — 0, 
we obtain, as expected, the two solutions Ei{k) = e{k) 
and E2 = Eq, describing respectively pure density and 
non interacting localized modes. A non-zero coupling 
hybridizes these two sets of excitations. Since X{k) and 
E{k) are also related by the relation (0) 



E'^{k)~ El = 2X{k)Eo 



(27) 



using the relations (|27| ) and (|2q ) , we obtain a two branch 
spectrum of dispersions 



E2 == 2Eo = 4A 



(28) 



These expressions are written in terms of e{k) and A, 
but they depend also exphcitely on X{k) through (H), 



which is a consequence of the dipolar approximation. As 
a result of the hybridization, the energy spectrum Ei(k) 
of the density fluctuations is shifted by a factor 2 towards 
lower energies relatively to the Feynman ansatz. There is 
in addition a new branch E2 of localized excitations at a 
constant energy which is twice the value of the bare local 
mode. It is of some interest to notice that it is possible 



to express the factors 



C(k) 
A(k) 



and 
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B(k) 



appearing in the 



coherent ground-state wavefunctions of the two branches 
(pq) in terms of the energies, namely 



C{k) 
A{k) 
D{k) 
B{k) 
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e{k) - E{k) 



e{k) + E{k) 
Eo - Eik) 



E„ + Eik) 



(29) 



For the special value k = krot which corresponds to 
the roton minimum, they are both equal, |C/A| (krot) — 
\D/B\ (krot) — 1/3, a consequence of the relation ([2q). 
The roton excitation then appears as involving the max- 
imum mixing of localized and density modes in equal 
proportions. 

We want now_ tn_compare these results with the ex- 
perimental dataa'tlEj obtained for the energy , spe ctoim 
E{k) (fig. 6) and for the structure factor ^(fc)^^^'^^ at 
two different pressures (fig. 7). Around the roton mini- 
mum, and over quite a large range of momentum, the 
relation ( |28| ) is in good agreement with the experimen- 
tal results and gives a spectrum which is reduced from 
the Feynman result by a factor of nearly two. It is 
interesting to recall that Feynman himself in his origi- 
nal paperEI noticed the factor 2 discrepancy between his 
ansatz and the experimental results around the roton 
minimum. Here we deduce it from a closed analytical 
solution. This interpretation of the roton excitation as 
resulting from the dipolar hybridization of two separate 
excitations ig-to be compared with the approach of Glyde 
and Griffin E3i3 which uses a dielectric formalism to de- 
scribe hybridized phonons and free-particle excitations. 
The main discrepancy is obtained in the low momentum 
region (fc < lA ), i.e. below the maxon momentum. 
As fc — > the experimental data is close to the Feyn- 
man result i.e. twice the value we obtain in (|28|). This 
regime is beyond the dipolar approximation we consider 
and sets the limit of validity of this approach. To be able 
to describe the linear phonon regime from the present 
approach, we should include multipolar terms. 



VI. HIGH ENERGY EXCITATION BRANCH 

A. Vortex loops 

The relation (£§) gives a second branch of excitations 
at the constant energy E2 = 4A. It describes a local- 
ized (dispersionless) excitation of energy which is twice 
the bare vortex core (local mode) energy. The intuitive 



physical picture of this mode is therefore of an excita- 
tion made out of two vortex-core excitations. Such an 
excitation in three dimensions may be viewed as a mi- 
croscopic vortex-loop, or a localized defect. The radius 
i? of a vortex-loop can be estimated using a Feynman- 
type forpmla for the energy of the associated circulating 
current EZl 



o 2 ^7 I "" 

^vortex = 27r p r-/n | — 



(30) 



where at T = 0, we take the density of the superfluid 
Ps to be the total density p and a is tlaa core size equal 
to the atomic radius namely a ~ 1.4v4E3. The radius R 
obtained from (BG) and which corresponds to E^o^tex — 
4A = 34.4if is R~ 5.ll. 

Experimental support for the interpretation of the hy- 
bridized localized state at E = 4A as an intrinsic exci- 
tation |Of the fluid is provided by critical velocity exper- 
imentsES. In phase-slippage studies of the critical veloc- 
ity through an orifice, the critical velocity is driven by 
the thermal nucleation of vortex-loops. The correspond- 
ing activation energy Ey is determined by the nucleation 

rate F given by the Arrhenius law F = Toexp I -^-^ ) and 

turns out to hec3Ey ~ 33±5K a value indeed very close 
to 4A = 34:.4:K. Moreover, the upper critical velocity Vc 
may Ipe, estimated to be the velocity of the vortex-loop 
e. Vr, = T;^ln f-^) ~ 20toZ^ a value close to 
ocity Ea. 
on cavitation in super- 



itselflji.e. v^ = ^ln(^ 

^ 2mR \ a / 

the largest measured critical 



Finally, recent experimentsEj 
fluid *He are analyzed in terms of thermal a nucleation 
process. The activation energy found is of the order of 
4A, indicating a possible connection between a micro- 
scopic vortex-loop and the thermal nucleation of cavita- 
tion. 

We also point out that this interpretation of the 
vortex-loop branch as distinct from the phonon-roton 
branch is further supported experimentally by ion ve- 
locity measureipaentalj and theoretically using hydrody- 
namic methodsca. 



B. High energy neutron scattering 

Recent high resolution neutron scattering data show 
that at high energy and momentum transfer, compared 
with the phonon-roton branch, the position of the scat- 
tering peakcj is consistently at a lower energy than the 
free-recoil value. At_such high energies the atoms of the 
liquid behave freelyO, therefore with a quadratic disper- 
sion Efree{k) — fi k'^ /2m. There is a marked broadening 
of the width of this peak where the free-recoil spectrum 
crosses the second branch of energy 4A at momentum 
kc — 2.iA~^. This data can be understood as resulting 
from the superposition of the free-recoil spectrum and 
the process of vortex-loops creation. Such a possibility 
for explaining the increased width of the dispersion curve 



using vortex-loop creation was suggested by Cowley and 
Woodalj. A recoiling atom can create vortex-loops and 
loose the corresponding energy and momentum, so that 
its remaining energy is 



E^{k) = 4nA- 



h^ {k - kc 
2to 



(31) 



where n is the number of emitted vortex-loops and kc^n 
is given by the equality of the free-recoil energy and the 
?ith vortex-loop energy, so that /cc,n — 2.3^/nA~^. We 
show the dispersion for the first two emitted vortex-loops 
in fig.8. Then, by substracting frorruthe experimental 
neutron-scattering data a LorentzianEl centered around 
the free-recoil energy and fitted to the highest energy 
tail, we find an additional scattering peak around the 
energy given by (pl|). In fig. 9 the experimental resultsBS 
are plotted, which correspond to fc's such that n = 1 or 2. 
At higher energy transfer the free-recoil spectrum can be 
split again at twice 4A and so on, as is indicated in fig.9d. 
There is a series of such splittings each time the recoiling 
atom has an energy which is aiujjiteger multiple of 4A. In 
fig. 10, the experimental resultscj for the high momentum 
transfer k — 10A~^ are shown. At this momentum the 
free- recoil energy is more than 17 times the vortex-loop 
energy of 4 A, and we find that the extra scattering is 
well described as a 'band' of splitted levels below the 
free-recoil spectrum. 



VII. CONCLUSION 

In this paper, we have presented a phenomcnological 
strong coupling theory for the superfluid phase of ^iJe in 
order to describe features that correspond to excitations 
of energies beyond the phonon part of the spectrum. This 
approach is based on the assumption that at short length 
scales, there exists a set of interacting localized excita- 
tions. This is very much in the spirit of the assumption of 
two- level tunneling states in dielectric glasses. These ex- 
citations are described by an effective Hamiltonian which 
depends on the strength X{k) of the interaction between 
these modes. We used experimental results obtained in 
neutron scattering in order deduce the behaviour oiX{k). 
This approach gives a hand to calculate analytical expres- 
sions for various physical quantities like the scattering 
cross-section Z{k) for single particle excitations, and the 
shift in the ground state energy. In addition, it allowed 
us to derive a relation between the condensate fraction 
and characteristics of the energy spectrum namely the 
sound velocity and the roton energy A. This gives a way 
to measure directly the condensate fraction. Moreover, 
it establishes the critical nature of A and its relation to 
superfluidity in close connexion with recent experimen- 
tal results. Then, we derived an alternative description 
of the set of localized modes which allowed us to predict 
the existence of a new branch of excitations at energy 4A 
beyond the phonon-roton branch. This branch which has 



not yet been observed directly shows up indirectly and 
allows to interpret phase slippage exper iments and to get 
a hint for the puzzling broadening observed at high en- 
ergy and momentum transfer neutron scattering. Since 
we interpret this excitation as the smallest vortex-loop 
(or local defect) of the superfluid, our model offers a way 
of describing both the phonoupcoton and the localized 
vortex-loops within one schemeCj. 
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FIG. 1. The experimental energy spectrum at saturation 
vapor pressure [3,9,32] compared with the Feynman phonon 
spectrum (Eq.l) based on the experimental static structure 
factor S{k) [33-35]. Also marked is the roton energy A and 
the termination of the quasi-particle branch at energy 2A 
(circled point). 
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FIG. 2. The interaction matrix X{k) from Eq.(4), based 
on the experimental energy spectrum at (a) saturation vapor 
pressure [3,9,32] and (b) P=24 atm [3]. The bare energy of 
the local mode is taken to be Eo = 2A. Notice that X{k) 
tends to zero at the termination point of the spectrum at 
high momentum. 
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FIG. 3. Comparison between the experimental scatter- 
ing cross-section [3] ^(fc) of single quasi-particle excitations 
(points) at I.IK and the theoretical curve (solid line), ob- 
tained using the experimental data for the energy spectrum 
at saturated vapor pressure -E(fe) [3,9,32] in Eq.(12). 
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FIG. 4. The integrand U^ {Eik) - 2A)^ appearing in the 
calculation of the reduction in the ground state energy A_Bg 
in Eq.(18). (a),(b) are at saturation vapor pressure and 
P=24atm respectively. Notice the large contributions around 
the roton momentum. 
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FIG. 5. Comparison between the experimental condensate 
fraction at difTerent densities (solid circles) [23] and Eq.(15), 
using an effective mass for the local modes of me// = 2.3m 
(solid line). 
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FIG. 6. Comparison between the experimental energy 
spectrum [3,9,32] (points) and the energy E\{k) in Eq.(28) 
(solid line), where the structure factor S{k) is obtained from 
independent measurements [33-35] . (a) and (b) correspond re- 
spectively to the saturation vapor pressure and to P=24 atm. 
The dashed line at energy E2 = 4A indicates the position of 
the localized branch of excitations (vortex- loop). 
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FIG. 7. Comparison between the experimental structure 
factor S'(fc) [33-35] (solid circles) and expression Eq.(28) (solid 
line) for the same two pressures as in Fig. 6, where the energy 
Eik) is obtained from independent measurements [3,9,32]. 
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FIG. 8. The dispersion relations of the free recoil (n = 0) 
and the lowest two levels of recoiling atoms that correspond 
to vortex-loops emission (n = 1,2) (Eq.31). 
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FIG. 9. Comparison between the experimental scatter- 
ing profiles [29] (solid circles) and the free-recoil Lorentzian 
(heavy dashed-line) at different momenta: (a) k — 2.8^^^, 
(b) k = S.OA^S (c) k = 3.2l-\ (d) k = 4.0A"\ The 
remaining extra scattering is the heavy solid- line. The vor- 
tex-loop energy (4A) is indicated by the vertical dotted-line, 
the free-recoil energy by the vertical long-dashed-line, and 
the n = 1 splitted energy Eq.(31) by the vertical dashed-dot 
line. At relatively low momenta (a,b) there is still a noticable 
phonon-roton peak around the 2A energy. In (d) the heavy 
vertical lines indicate the n = 2 vortex-loop emission by the 
recoiling atom. 
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FIG. 10. (a) Experimental scattering profile [44] at 
k = 10. OA^^ (solid circles). The free-recoil Lorentzian (heavy 
dashes-line) and extra scattering (heavy solid-line) are plot- 
ted. The free-recoil energy is indicated by the long-dashed 
vertical line, (b) The spectrum of the free-recoil (dash-dot) 
and the 17 splittings due to vortex-loop excitations (Eq. 31) 
(alternating solid and dashed lines). 
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